Matrices coupled in a chain. II. Spacing functions 
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Abstract. For the eigenvalues of p complex hermitian n x n matrices coupled in a chain, 
we give a method of calculating the spacing functions. This is a generalization of the one 
matrix case which has been known for a long time. 



1. Introduction 

Let us recall here a few facts concerning the case of one matrix. For a n x n complex her- 
mitian matrix A with matrix elements probability density exp[— tr the probability 
density of its eigenvalues x := {x%,X2, x n } is [1] 



F(x) oc exp -J2 v ( 
3=1 

oc det[K(xi, Xj)]ij=i,... }n , 



l<i<j'<n 



l.la) 
(1.16) 



where V(x) is a real polynomial of even order, the coefficient of the highest power bein£ 
positive; K(x, y) is defined by 



K{x,y) := exp 



l -V(x) - lv(y) 



n—l 1 



'1.21 



J i=0 



Pi(x) is a real polynomial of degree i and the polynomials are chosen orthogonal with the 
weight exp[— V(x)], 

J Pi(x)Pj(x) exp[-V(x)}dx = hiSij. (1.3) 

Here and in what follows, all the integrals are taken from -co to +oo, unless explicitly 
stated otherwise. 
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The correlation function Rk{xi, Xk), i.e. the density of ordered sets of k eigenvalues 
within small intervals around x±, Xk, ignoring the other eigenvalues, is 



Rk(xi,...,Xk) ■= ^ J^, J F(x.)dx k +i...dx n 
= det[K(x i ,x j )] id=lj ... j k- 



(1.4) 



The spacing function E(k,I), i.e. the probability that a chosen domain I contains 
exactly k eigenvalues (0 < k < n), is 



E(k,I) 



n 



k\{n-k)\ 

k 



F(x) 



3 = 1 



n [j-xwi 



dxi...dx 7 



1 / d 



P(*,J) 



z=-l 



fc! V^, 

where x(x) is the characteristic function of the domain /, 

■v(x) ■= I lj if ^ G J ' 
v 10, otherwise, 



(1.5) 



(1.6) 



and R(z, I) is the generating function of the integrals over / of the correlation functions 
Rk(xi, ...,Xk), 



/n n 

j=l k=0 K - 

n, ^ k = o, 

Pk ~) J Rk(xi,...,Xk) Y\_x(xj)dxj, otherwise. 

I 3 = 1 

The R(z, I) of Eq. (1.7) can be expressed as a determinant 

R( z i -0 = det[Gy]j 5 j = o,..., n -i) 



(1.7) 



(1.8) 



(1.9) 



where, using the orthogonality, Eq. (1.3), of polynomials Pi (a;) and splitting the constant 
and linear terms in z, Gij reads 



G 



1 



t3 ~ K 



P l (x)P J (x) exp[-V(x)][l + z X (x)]dx = 5 tJ + G %3 . 
Finally, R(z, I) can also be written as the Fredholm determinant 

n 

R(zJ) = Y\[l + X k (zJ)] 



(1.10) 



(1.11) 



k=i 
2 



of the integral equation 



/ 



N(x,y)f(y)dy=Xf(x), 



(1.12) 



where remarkably the kernel N(x,y) is simply zK(x, y)x(y) with K(x,y) of Eq. (1.2). 
The Xi(z, I) are the eigenvalues of the above equation and also of the matrix [Gij]. 

These results can be extended to a chain of p complex hermitian n x n matrices. We 
consider the probability density for their elements 

tr j^i(^i) + V 2 {A 2 ) + ■■■ + K p _i(4 p _i) + \v p {A p ) 



JF(Ai, • • • , A p ) oc exp 

x exp [tr {ci A X A 2 + c 2 A 2 A 3 H h c v - X A v _ x A p }] 



(1.13) 



Here Vj(x) are real polynomials of even order with positive coefficients of their highest 
powers and the Cj are real constants such that all the integrals which follow converge. For 
each j the eigenvalues of the matrix Aj are real and will be denoted by xj := {xji, Xj 2 , 
Xj n }. The probability density for the eigenvalues of all the p matrices resulting from Eq. 
(1.13) is [2-5] 



F(xi; ...;x p ) = Cexp 



71 r 1 1 V 

S I 2 Vl<yXl ^ + V2<yX2r ^ H h V p-l( X p-lr) + ^V p (x pr ) I 

r=l ^ ' 



X 



W ("^Is x lr)i, x ps %pr) 
l<r<s<n 



det \e ClXlrX2 °] det [e C2X2 ' 



X3 S 



'det [g c p — 1 lra: 'P s j 

(1.14) 



= C 
where 



ps 



X 



pr ) 



Kr<s<n 



p-1 



"f[ det [w k (x kr , x fc +i s )] rjS=1) 



.k=i 



w k {x,y) := exp 



~\ v k{x) - ^V k+1 (y) +c k xy 



(1.15) 



(1.16) 



and C is a normalisation constant such that the integral of F over all the np variables Xi r 
is one. 

The correlation function 



Rk\,...,k p \ x lli •••'■> %pli •••■> x pkp) 



r p 

/ F(x i; ...;x p ) 
J i=i 



n! 



V V r J= k j+ l 



(1.17) 



was calculated in a previous paper [6] to be an m x m determinant (m = k\ + ... + /c p ) 

-Rfci ,...,/cp (2-11? •••) •••> x pkp ) 

= det [Ky (X ir , X^)]^. 1 p;r 1 /.-.:, I • ( L18 ) 



This is the density of ordered sets of kj eigenvalues of Aj within small intervals around 
x jU x jkj for j = 1, 2, p. 

Here we will consider the spacing function E(ki, I\; ...; k p , I p ), i.e. the probability that 
the domain Ij contains exactly kj eigenvalues of the matrix Aj for j = 1, ...,p, < kj < n. 
The domains Ij may have overlaps. As in the one matrix case one has evidently 



1 f d 
E{k u h- : ...- : k p ,Ip) = — 



1 / d 



kp\ \ 



R(zi,Ii', z p , I p ) 



k\\ \dzi 
with the generating function 

/V n 
F(xi;...;x p ) JJ JJ [1 + Zj Xj (x jr )] 
A 1 1 



dx 



j f 



= E-E 



j=l r =l 

Pk!,...,k p fci ^fc p 

fci!...ife p ! 1 " p ' 



(1.19) 

(1.20) 
(1.21) 




I j dx j f 



3 r=l 



— ••• — kp — 0, 
Rk 1 ,...,k p (xn, • • • , x lkl ; x p i, • • • , x p k v ), otherwise, 

(1.22) 



Xj( x ) being the characteristic function of the domain Ij, Eq. (1.6). 

The function R(zi, I\ \ ...; z p , I p ) will be expressed as a n x n determinant. It will also 
be written as a Fredholm determinant, the kernel of which will now depend on the variables 
zi, z p and the domains I±, I p in a more involved way than in the one matrix case. 
In particular, it does not have the remarkable form mentioned after Eq. (1.12). 



2. The generating function R(z\, I\; ...; z p , I p ) 

The expression of F, Eq. (1.15), contains a product of determinants. As the product of 
differences 

A(xi)= Y[ (x ls -x lr ) (2.1) 

l<r<s<n 

and det[w\(xi r ,X2 S )] are completely antisymmetric and other factors in the integrand 
of Eq. (1-20) are completely symmetric in the variables in, xi n , one can replace 
det[wi (xi r , X2 S )] under the integral sign in Eq. (1.20) by a single term, say the diagonal one, 
and multiply by n\. This single term is invariant under a permutation of the variables x\ r 
and simultaneously the same permutation on the variables xir- Therefore, after integration 
over the the integrand, excluding the factor det[w2(^2r, X3s)], is completely 
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antisymmetric in the variables X21, xm and so one can replace the second determinant 
det[w2(x2r, %3s)] by a single term, say the diagonal one, and multiply the result by n\. In 
this way, under the integral sign one can replace successively each of the p — 1 determinants 
det[wk(xkr , Xk+is)] by a single term multiplying the result each time by n\ 



R(z u h;...;z p J p ) = (nlf-'C J A( Xl )A(x p ) 

V n 

n n [ i+z jXj(xjr) 



p—l n 

YlYl w j( x jr,Xj+lr) 

3 = 1 r=l 



3 = 1 r =l 



dx 



3 r 



(2.2) 



Recall that a polynomial is called monic when the coefficient of the highest power is 
one. Also recall that the product of differences A(x) can be written as a determinant 



A(x) = det[x{ x ] = det^-ifo)] = det[g,_!(^)], 



(2.3) 



where Pj(x) and Qj(x) are arbitrary monic polynomials of degree j. As usual, we will 
choose these polynomials real and bi-orthogonal [6] 



J P J (x)(w 1 *w 2 * ... * w p -i)(x,y)Q k (y)dxdy = hj8 jk , 



(2.4) 



with the obvious notation 



(f*9)(x,y)= / f(x,H)g(t;,y)dC 



(2.5) 



The normalization constant C is [6], 

C=(nl)-*I[hr 1 . 



n-l 



(2.6) 



i=0 



Now expand the determinant as a sum over n\ permutations (i) := ( °' ' ) , ir(i) 



being its sign, 



j 7 



det[P s _i(xi r )] = J^7r(i)P il (xii)P i3 (a;i2)...Pi„(xi n ). 



(2.7) 



(0 



Doing the same thing for det[Q s -i(x pr )] and integrating over all the np variables xj r ; 
j = 1, ...,p; r = 1, n in Eq. (2.2), one gets 



R(z 1 ,I 1 ;...;z p ,I p ) = ^^2^Tv(i)ir(j)G njl G i2j2 ...G lnJ 

(0 Cj) 



(2.8) 



where 



■p-i 

j J W k {x k ,Xk+l) 

.fc=l 



Qj(x p ) 



Y[[l + ZkXk(xk)]dx k 



.fc=i 



(2.9) 



When all the Zfc vanish, is equal to 5^ as a consequence of the bi-orthogonality, 
Eq. (2.4), of the polynomials Pi(x) and Qi(x). Let us define Gij as follows 



G 



*3 



Gij 5ij , 



so that 



1 



Gij = ^ / P*(xi) 



p-i 

j J w k (x k ,x k+1 ) 



Qj(x p ) 



Y[[l + z k Xk(xk)} ~ 1 



k=l 



{[ dxy 



k=l 



(2.10) 



(2.11) 



Any n x n determinant is the product of its n eigenvalues and therefore one has 

n 

R(zi, h] z p , I p ) = J J [1 + Afc(zi, ii; ...; z p , ip)] , 



(2.12) 



k=i 



where the Xk(zi, h; ...; z p , I p ) are the n roots (not necessarily distinct, either real or pairwise 
complex conjugates, since is real) of the algebraic equation in A 



det[Gij - XSij] = 0. 



(2.13) 



One can always write a Fredholm integral equation with a separable kernel whose eigen- 
values are identical to these (cf. reference [7] for the case of p = 2 matrices). Indeed, for 
any eigenvalue A the system of linear equations 



n-l 

Gijij = Ac 

3=0 



(2.14) 



has at least one solution i = 0, n — 1, not all zero. Multiplying both sides of the above 
equation by Qi(x), summing over i and using Eq. (2.11) gives the Fredholm equation 



N(x,Xp)f(x p )dx p = Xf(x), 



(2.15) 



where 
f(x) 

N(x,x p ) 



n-l 



i=0 

^li Qi(X) j Pt(Xl) 



(2.16) 



p-i 

Y[ Wk(x k ,X k +l) 
.k=l 



Yli 1 + z kXk(x k )] - 1 



k=l 



■p-1 



n dxk 



.fc=i 



(2.17) 
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Hence if A is an eigenvalue of the matrix [Gij], it is also an eigenvalue of the integral 
equation (2.15). Conversely, since the kernel N(x,x p ) is a sum of separable ones and 
since the polynomials Qi(x) for i = 0, n — 1 are linearly independent, if A and f(x) 
are, respectively an eigenvalue and an eigenfunction of this integral equation, then f(x) is 
necessarily of the form 

n-l 

= (2.i8) 

i=0 

and the i = 0, n — 1, not all zero, satisfy Eq. (2.14). Therefore A is a root of Eq. 
(2.13). 

When one considers the eigenvalues of a single matrix anywhere in the chain, disre- 
garding those of the other matrices, everything works as if one is dealing with the one 
matrix case and formulas (1.2), (1.5), (1.7) and (1.11) are valid with minor replacements. 
Similarly, when one considers properties of the eigenvalues of k (1 < k < p) matrices 
situated anywhere in the chain, not necessarily consecutive, everything works as if one is 
dealing with a chain of only k matrices; the presence of other matrices modifying only the 
couplings. 

To say something more about the general case seems difficult. 

When Vj(x) = cijX 2 , j = 1, p, then the polynomials Pj(x) and Qj(x) are Hermite 
polynomials Pj(x) = Hj(ax), Qj(x) = Hj(f3x), the constants a and (3 depending on the 
parameters cij and the couplings Cj. In this particular case the calculation can be pushed 
to the end (see appendix). 



Appendix 

For Vj(x) = a,jX 2 : j = 1, p, setting 



11 a j,. c {,r. y) := exp ( -^ax 2 - h)y 2 + cxy j . (-4.1 ) 



one gets according to Eq. (2.5) the multiplication law 

1/2 



(W aAc * W a / )6 / )C /)(x,y) = ( — —J W a ",b",c"(x,y), (A.2) 
where 

c 2 c' 2 cc' 

b + a' b + a' b + a' 

From Eq. (1.16), Wk(x,y) = W akjak+ljCk (x,y) and a repeated use of the above multiplica- 
tion law yields 

W(x, y):=(w 1 *w 2 *---* u> p _i)(x, y) = d x W a ^ c (x, y), (AA) 
where a, 6, c and d are constants depending on the parameters ai, a p and Ci, Cp—±. 
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The orthogonality relation (2.4) of the polynomials Pj(x) and Qj(x) takes the form 

J Pj(x)W(x,y)Q k (y)dxdy = hj5 jk , (A.5) 

namely the same relation as in the two matrix case with the weight W(x, y), an exponential 
of a quadratic form in x and y. It follows that Pj(x) and Qj(x) are Hermite polynomials 
of x times a constant 

Pj(x) = Hj(ax), a : = (^^) ' \ (A.6) 

/ab-c 2 \ 1/2 

Q j (x) = H j (px), f 3 --={-^) ; ( A7 ) 

hj = 7Zi (4=y^^-d. (A8) 



(a6 _ c 2)i/2 

The eigenvalue density of the matrix Ai, for example, ignoring the eigenvalues of other 
matrices, is from Eq. (1.18) 

n-1 



1 f 

R^x) = K n (x, x) = ^2 T-Pj{x) / W{x,y)Q j {y)dy 

3=0 3 J 



3'- 

which in the large n limit is a semi-circle of radius \/2n/a. Thus in this particular case 
of coupled matrices one recovers Wigner's "semi-circle law" for the eigenvalues of a single 
matrix. 

The kernel of the integral equation (2.15) is given by Eq. (2.17) with Pj(x), Qj(x) and 
hj as given above. To get further, one has to take explicitly the domains 7, into account. 
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